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Abstract 



It is noted that two-by-two "S" matrices in multilayer optics can be repre- 
sented by the Sp(2) group whose algebraic property is the same as the group 
of Lorentz transformations applicable to two space-like and one time-like di- 
mensions. It is noted also that Wigner's little groups have a sliderule-like 
property which allows us to perform multiplications by additions. It is shown 
that these two mathematical properties lead to a cyclic representation of the 
S-matrix for multilayer optics, as in the case of ABCD matrices for laser cav- 
ities. It is therefore possible to write the iV-layer S-matrix as a multiplication 
of the N single-layer S-matrices resulting in the same mathematical expres- 
sion with one of the parameters multiplied by N. In addition, it is noted, as in 
the case of lens optics, multilayer optics can serve as an analogue computer for 
the contraction of Wigner's little groups for internal space-time symmetries 
of relativistic particles. 
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I. INTRODUCTION 



In our previous paper on multilayer optics [1], it was shown that the complex two-by-two 
S-matrix formalism is equivalent to a two-by-two real matrix representation of the Sp(2) 
group, which shares the same algebraic property as the Lorentz group applicable to two 
space-like and one time-like dimensions. This group has three independent parameters. 
It was shown furthermore that, under certain conditions, one of the off-diagonal elements 
vanishes, and the three remaining elements can be computed analytically. We called this 
the Iwasawa effect [1]. In this paper, we remove those "certain conditions" and achieve the 
same kind of simplification for all possible multilayer cases. 

Indeed, the group Sp(2) plays the central role in both quantum and classical optics, 
including multilayer optics [2] . It consists of two-by-two real matrices whose determinant is 
one. Each matrix contains at most three independent parameters. It is thus a simple matter 
to multiply two or three matrices. However, for multiplication of a large number of matrices 
presents a new problem. The product of those many matrices will also be one two-by-two 
matrix with a unit determinant, but how can we calculate their elements? 

For example, let us look at laser cavities. It consists of a chain of N identical two-lens 
systems, where N is the number of cycles the light beam performs. The resulting ABCD 
matrix can be written as a multiplication of N identical matrices, but the resulting matrix 
has the same mathematical form as that for the single cycle [3]. 

Can we then expect a similar cyclic property in multilayer optics? We have shown in 
Ref. [1] that the N-dependence can be made quite transparent if the multilayer S-matrix [4] 
is reduced to the Iwasawa form. In this paper, we present the cyclic property for the most 
general form of multilayers, without the restriction we imposed in our previous paper [1]. 
We shall show that the core of the S-matrix takes the form 

cos a — sincA /cosh/9 sinh/9\ /l 0\ (1 j\ , . 

sina cosaj' \smhf3 cosh/^J ' U \) ' ° T \0 l) ' () 

These matrices form the core of Wigner's little groups applicable to the internal space-time 
symmetries of relativistic particles [5,6]. We note here that these matrices have the following 
interesting property. 

We cannot write (cosai x cosa 2 ) = cos(o;i + a 2 ) because it is wrong. However, in the 
two-by-two matrix form, 



cos«i — sinaA /cosa 2 _ sina: 2 \ _ ( cos(o;i + a 2 ) — sin(a:i + a 2 ) 
sinai cosai / \sina: 2 coso; 2 / \ sin(«i + a 2 ) cos(a;i + a 2 ) 



(2) 



and similar expressions for the remaining matrices in Eq.(l). We call this the sliderule 
property of Wigner's little groups. 

If they are cycled N times, they take the form 

cos(iVa) -sin(iVa)\ /cosh(7V/3) sinh(iV/3) \ 
sin(iVa) cos(iVa) ) ' I sinh(/3) (cosh/3) ) ' 



1 0\ j f \ 

and 



n 1 ij> """" \o i y 
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respectively. This mathematical instrumentation works for laser cavity optics [3]. The 
question is whether this is applicable to multilayer optics. 

The purpose of this paper is to show that the answer to the above question is YES. 
We note first that the S-matrix consists of iV cycles. Each cycle consists two phase-shift 
matrices, one boundary matrix and its inverse, and this cycle does not take any of the forms 
given in Eq.(l) if we start the cycle from the boundary. In this paper, we show that it is 
possible to obtain the core in the form of Eq.(l) if we start the cycle from somewhere within 
one of the media between the two boundaries. 

Throughout this paper, we avoid group theoretical languages and rely on explicit two-by- 
two matrices with real elements. However, in so doing, we are going through an important 
group theoretical aspect which became known to us only recently, namely on contractions 
of Wigner's little groups. This aspect was discussed in detail in a recent paper on lens 
optics [7]. Thus, we shall borrow some of the mathematical identities from that paper. 

In addition, in the present paper, we observe that Wigner's little group has sliderule 
properties which allow us to convert multiplications into additions. This property was noted 
for one of the little groups in the paper of Han et al. In this paper, we shall show that all 
three of the little groups have the same sliderule property, using Eq.(2). 

In Sec. II, we formulate the problem in terms of the S-matrix method widely used in 
multilayer optics [4,8,9], and show that the complex S-matrices can be transformed to real 
matrices by a conjugate transformation, and thus to the algebra of the Sp(2) group which 
is by-now a familiar mathematical language in optics. In Sec. Ill, we import from the 
literature mathematical identities useful for the purpose of the present paper. They are 
derivable from Wigner's little groups and their contractions. In Sec. IV, using the cyclic 
property of Eq.(3), it is shown possible to write the multilayer S-matrix as a multiplication 
of the N single-layer S-matrices resulting in the same mathematical expression with one of 
the parameters multiplied by N. In Sec. V, it is pointed out that the mathematical identities 
presented in this paper can be tested experimentally. We discuss the condition under which 
the system can achieve the Iwasawa effect [1]. 



II. FORMULATION OF THE PROBLEM 

It was noted our previous paper that one cycle in iV-layer optics starts with the boundary 
matrix of the form [10] 

T3( n \ - ( cosh(?7/2) sinh(?7/2) \ 
W Wnh(77/2) cosh(??/2) J ' { ) 

which describes the transition from medium 2 to medium 1, taking into account both the 
transmission and reflection of the beam. As the beam goes through the medium 1, the beam 
undergoes the phase shift represented by the matrix 

When the wave hits the surface of the second medium, the corresponding matrix is 
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cos W 2 ) -sinh(r//2)\ ( . 

1 V) v -sinh(V2) cosh(?7/2) )' W 

which is the inverse of the matrix given in Eq.(4). Within the second medium, we write the 
phase-shift matrix as 

P(<h) = ( e Q J 2/ 2)- (7) 

Then, when the wave hits the first medium from the second, we have to go back to Eq.(4). 
Thus, one cycle consists of 



Mi 



cosh(r?/2) sinh(?7/2) \ /e-^ 1 / 2 \ / cosh(^/2) -sinh(??/2)\ 
sinh(r]/2) cosh(r?/2) J \ e^ 1 / 2 J \ - sinh(?7/2) cosh(r7/2)J 

x{ o e^J- (8) 

This matrix contains complex numbers, but we are interested in carrying out calculations 
with real matrices. This can be done if we make the following conjugate transformation [1] 
Let us next consider the matrix 

Then we have shown in our previous paper that 

M 2 = C Mi C~\ (10) 



with 



o Wcos(^/2) -sin(&/2)\ /e""/ 2 
• 7 ' ? _l e-^ 2 Asin(0i/2) cos^/2) J V e^ 2 

/cos(0 2 /2) -sin(0 2 /2)\ (n) 



sin(0 2 /2) cos(0 2 /2) 

The conjugate transformation of Eq.(10) changes the boundary matrix B(r)) of Eq.(4) 
to a squeeze matrix 

^=(T e-/ 2 )' (12) 

and the phase-shift matrices P(<j>i) of Eq.(5) and Eq.(7) to rotation matrices 

r>/A \ _ fcos(&/2) -sin(0i/2)\ 

Uin(^/2) cos(^/2) J' ^ 

with i — 1, 2. 

Indeed, the matrices Mi and M 2 can be written as 
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M l = B{ V )P{<p 1 )B(-i 1 )P{<t> 2 ), 



M 2 = S(r ] )R( ( j )l )S(-r ] )R( ( f )2 ). (14) 

The matrix M 2 can be obtained from Mi by the conjugate transformation in Eq.(lO). Con- 
versely, Mi can be obtained from M 2 through the inverse conjugate transformation: 

Mi = C- 1 M 2 C. (15) 

In addition, the conjugate transformations have the following properties. 

(M 2 f = C (M.f C-\ (M 1 f = C' 1 (M 2 f C. (16) 

Thus, we can study M 2 in order to study Mi. The advantage of M 2 is that it consists of 
real matrices. The group of these matrices is called Sp(2) which is like (isomorphic) the 
Lorentz group applicable to three space and one time dimensions. This group contains very 
rich group theoretical contents including those of Wigner's little groups. We intend to study 
M 2 in terms of those little groups. 

The problem is that M 2 takes a simple form, and (M 2 ) 2 is manageable, but we cannot 
predict what form {M 2 ) N takes. In this paper, we shall construct the core matrix of the 
form of Eq.(l) for multilayer optics. Then, as we can see in Eq.(3), the chain effect is 
straight-forward. We shall calculate (M 2 ) N first and then (Mi)^. 



III. MATHEMATICAL IDENTITIES FROM THE LORENTZ GROUP 

Wigner's little groups were formulated for internal space-time symmetries of relativis- 
tic particles [5,6]. However, they produced many mathematical identities useful in other 
branches of physics, including classical layer optics which depends heavily on two-by-two 
matrices. The correspondence between the two-by-two and four-by-four representations of 
the Lorentz group has been repeatedly discussed in the literature [1,3,7]. In the two-by-two 
representation, we write the rotation matrix around the y axis as 



cos(0/2) -sin(0/2)\ 
, sin(0/2) cos(0/2) J ' 

and the boost matrices along the z and x axes as 

e r >l 2 \ / cosh A sinh A 

e- r '/ 2 J ' I sinh A cosh A 



(17) 



(18) 



respectively. We shall use only these three matrices in this paper. 

We use the following identity which Baskal and Kim introduced recently in their paper 
on lens optics and group contractions [7,11]. 

' e^/ 2 \ / cos (0/2) -sin(0/2)\ (e-' r >/ 2 
e-^l 2 ) \ sin(0/2) cos(0/2) J \ e^ 2 



cos(#/2) -sin(0/2)\ /cosh A sinh A \ /cos(0/2) -sin(0/2)\ 
sin(0/2) cos(#/2) J \ sinh A cosh A J Un(0/2) cos(#/2) ) 



(19) 
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with 



cos(0/2) = cosh A cos 9, 

r>„ cosh A sin 9 + sinh A 

e v = . (20) 

cosh A sin 9 — sinh A 

The left-hand side of the above expression is one rotation matrix sandwiched by one boost 
matrix and its inverse, while the right-hand side consists of one boost matrix sandwiched 
between two identical rotation matrices. 

The left-hand side of Eq.(19) is the same as the first three matrices of the core matrix 
M 2 given in Eq.(ll). However, the fourth matrix is a rotation matrix. Since one-rotation 
matrix multiplied by another rotation matrix is still a rotation matrix, the core matrix M 2 
is one boost matrix sandwiched between two different rotation matrices. Thus, the problem 
is to find a transformation which will make those two rotation matrices the same, and go 
back to the form of the left-hand side of Eq.(19). We shall come back to this problem in 
Sec. IV. 

If we complete the matrix multiplications of both side, the result is 



cos(0/2) — e v sin(0/2) \ _ / cosh A cos 9 —(cosh A sin 9 + sinh A) 
e _?? sin(0/2) cos(0/2) / \ cosh A sin 9 — sinh A cosh A cos 9 



(21) 



Then we can write <fi and rj in terms of A and 9 as given in Eq.(20). The parameters A and 
9 can be written in terms of and r] as 

cosh A = (cosh 77)^1 — cos 2 (0/2) tanh 2 ??, 

cos 9 = C ° SW2) (22) 

(cosh ??) yl — cos 2 (0/2) tanh 2 rj 

The above relation is valid only for (cosh A sin 9/2 — sinh A) is positive. If it is negative, 
the left-hand side of the above expression should be 

e^ 2 \ / cosh(x/2) - sinh(x/2) \ / e'^ 2 
e-^ 2 ) V - sinh(x/2) cosh(x/2) ) { e^ 2 



cosh(x/2) -e T 'sinh(x/2)\ 
-e^ sinh(x/2) cosh(x/2) J' 

with 

cosh(x/2) = cosh A cos 9, 



(23) 



sinh A — cosh A sin 9 
Conversely, A and 9 can be written in terms of x an d f] as 



9 „ cosh A sin 9 + sinh A 
e = : — . (24) 
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cosh A = (cosh r?) y cosh 2 (x/2) — tanh 2 ?], 



a cosh(x/2) 

cost/ = . . (25) 

(cosh 77) \J cosh 2 (x/2) — tanh 2 ^ 

An interesting case is when sinh A — cosh A sin 9 becomes zero, and 77 becomes very large. 
If we insist that 

e" sin(0/2) = u, (26) 
remain finite, then 0/2 must become very small. On the right-hand side, 

u — 2 sinh A, wit/i sin = tanh A. (27) 
The net result is that both sides take the form 

I - 2 f hA ). (28) 

In their recent paper [7], Kim and Baskal studied in detail the transition from Eq.(21) to 
Eq.(23) through Eq.(28), and showed that the one- lens camera goes through this transition 
as we try to focus the image. Mathematically, the system goes through group contraction 
processes. In the present paper, we show that the same contraction process can be achieved 
in multilayer optics. 

IV. CYCLIC REPRESENTATION OF THE S MATRIX 

It was noted in Sec. II that each cycle consists of 

(SR.S- 1 ^) , (29) 

with 

R x = R{<f> x ), R 2 = tf(0 2 ), (30) 

of Eq.(13) respectively. The squeeze matrix S is given in Eq.(12). For the layer consisting 
of N cycles, let us consider the chain 

= ( y SR 1 S- 1 R 2 ) (SR 1 S- 1 R 2 ) ( y SR 1 S- 1 R 2 ) .... (SR^Ri) . (31) 

According to Eq.(19), we can now write SRiS' 1 in the above expression as 

SR^- 1 = R 3 X R 3 , (32) 

with 

_/cos(0 3 /2) -sin(0 3 /2)\ /cosh A sinh A \ , . 

3 Wn(0 3 /2) cos(0 3 /2) )' UinhA cosh A/' 1 J 



and 

cosh A = (cosh r])\Jl — cos 2 (0i/2) tanh 2 77, 

cos(0i/2) , . 

cos0 3 = , = ' ==• (34) 

(cosh 77)^1 — cos 2 (0i/2) tanh rj 

The parameters A and <p 3 are determined from r\ and 0i which are the input parameters 
from the optical properties of the media. 
The chain of Eq.(31) becomes 

M 2 N = (R 3 XR 3 R 2 ) (R 3 XR 3 R 2 ) (R 3 XR 3 R 2 ) .... (R 3 XR 3 R 2 ) . (35) 

Let us next introduce the rotation matrix R(a) as 

R(a) = (R 2 ) 1/2 R 3 , (36) 

with 

OL = 03 + ^02, (37) 

where 2 is an input parameter. Since 3 is determined by r\ and 0i, the rotation angle a is 
determined by the three input parameters, namely i], 0i, and <p 2 . 
In terms of R = R(a), the chain of Eq.(35) becomes 

Mf = R 3 R~ 1 (RXR)(RXR)(RXR)....(RXR)R- 1 R 3 R 2 . (38) 

Since i^iT 1 = R 2 1/2 and R^R 3 R 2 = R l 2 ' 2 from Eq.(37), 

Mf = (R 2 )~ 1/2 (RXR)(RXR)(RXR) (RX R)(R 2 ) 1/2 . (39) 

According to Eq.(19) and Eq.(21), we can now write RXR as 

j^Xft _ ( cosh ^ cos a ~ (cosh A sin a + sinh A) \ ,^ . 

V cosh A sin a — sinh A cosh A cos a ) ' 

According to the formulas given in Sec. Ill, especially Eq.(19), RXR can also be written as 

RXR = ZAZ~\ (41) 

with 

e-^J- ^ 

Now the two-by-two matrix A can take one of the following forms. 

If the off-diagonal elements of the matrix of Eq.(40) has opposite signs, the A matrix 
becomes 
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/cos(0/2) -sin(0/2)\ 

Um(0/2) cos(0/2) { 6) 

with 

cos(0/2) = cosh A cos a, 

9 t cosh A sin a + sinh a , , 

e ? = . (44) 

cosh A sin a — sinh A 

If, on the other hand, the off-diagonal elements of the matrix RXR have the same sign, the 
matrix A should be written as 



_ / cosh( X /2) -sinh( X /2)\ 

V-sinh( X /2) cosh( X /2) )> ^ 

with 

cosh(x/2) = cosh A cos a, 

2 £ cosh A sin a + sinh A 
sinh A — cosh A sin a 

We note from Eq.(43) and Eq.(45) that the matrix A takes circular or hyperbolic forms 
depending on the sign of the lower- left element of Eq.(40) which is 

sinh A — (sin a) cosh A, (47) 

and note that this expression can become from a positive to negative number continuously 
as the parameters A and a vary. These two parameters are determined from the reflection 
and transmission properties of the media. 

While expression of Eq.(47) makes the continuous transition, it has to go through zero. 
If it vanishes, 

RXR=(l " 2s j nhA ). (48) 

The transition of A from Eq.(43) to Eq.(45) through this process has been discussed in detail 
in Ref. [7] in connection with the contraction of Wigner's little groups. 

As we noted in Sec. II, the matrix A has the desired cyclic property. Thus, 

M? = (R 2 )~ 1/2 [(ZAZ- 1 ) (ZAZ- 1 ) (ZAZ- 1 ) ... (ZAZ' 1 )] (R 2 ) 1/2 . (49) 

Consequently 



M 2 N = (R 2 y 1/2 Z A N Z- 1 (R 2 ) irz . (50) 



1/2 



If A takes the form of Eq.(43), 



N _ /cos(iV0/2) -sin(N<j>/2)\ , , 

"~ \sm(N<j)/2) cos(iV0/2) J' {0 ) 
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For A given in Eq.(45), 



A 



N 



cosh(iVx/2) - sinh(iVx/2) \ 
-sinh(A^x/2) cosh(iVx/2) ) 



As Eq.(48), 



(RXR) 



N 



1 -2N sinh A 
1 



(52) 



(53) 



Then the calculation of (M 2 ) N for the iV-layer case is straight-forward. We can now 



compute the matrix {M\) N using the conjugate transformation of Eq.(16). Let us write our 
result in two-by-two matrices: 



M 2 N = 



x 



cos(0 2 /4) - sin(0 2 /4) 
sin(0 2 /4) cos(0 2 /4) 



e^ 2 
e~«/ 2 



cos(iV0/2) - sin(iV0/2) \ 
sin(A^>/2) cos(A^>/2) J 



o 



e«/ 2 



cos(0 2 /4) sin(0 2 /4) 
- sin(0 2 /4) cos(0 2 /4) 



(54) 



for A of Eq.(43). For A of Eq.(45), 

•cos(0 2 /4) -sin(0 2 /4)\ /e«/ 2 
sin(0 2 /4) cos(0 2 /4) J { e^ 2 



cosh(A^x/2) - sinh(iVx/2) \ 
- sin(A^x/2) cos(iVx/2) ) 



x 



o 



cos(0 2 /4) sin(0 2 /4) 
sin(0 2 /4) cos(0 2 /4) 



(55) 



If the lower-left element given in Eq.(47) vanishes, we have to go back to Eq.(39) and 
Eq.(48), and write 



M 2 N 



cos(0 2 /4) - sin(0 2 /4) \ / 1 -2iVsinhA 
sin(0 2 /4) cos(0 2 /4) )\0 1 



cos(0 2 /4) sin(0 2 /4) 
- sin(0 2 /4) cos(0 2 /4) 



(56) 



As we noted in Sec. II, we use M 2 and M 2 for mathematical convenience. In the real 
world, we have to use Mi and . It is not difficult to write this expression using the 
conjugate transformation of Eq.(15). It can be written as 



N 



-#2/4 







,102 /4 



cosh(£/2) sinh(£/2) \ 
sinh(£/2) cosh(£/2) ) 



-iN</>/2 







iN<f>/2 



X 



cosh(£/2) -sinh(£/2)\ /e^ 2 / 4 
- sinh(£/2) cosh(£/2) ) { e~^ 2 / 4 



(57) 



if A takes the form of Eq.(43) with a positive value of Eq.(47). If it takes the form of Eq.(45) 
with a negative value of Eq.(47), 

• e -^ 2 / 4 \ /cosh(£/2) sinh(f/2)\] / cosh(A^x/2) i sinh(7Vx/2) \ 
e^ 2 / 4 J Uinh(^/2) cosh(£/2) sinh(A^/2) cosh(N X /2) ) 



N 



X 



cosh(£/2) - sinh(£/2) \ / e i<fe / 4 
-sinh(f/2) cosh(£/2) J { e"^ 2 / 4 



(58) 
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If the expression of Eq. (47) vanishes, 



M N_fe' i(p2/4 \fl-iNwahX iiVsinhA \ /e^ a/4 \ . ■ 

1 "I e^/yK -iiVsinhA 1 + iTVsinhAjv e^/ 4 J ' ^ 

This is not yet the S-matrix. The first and the last layers have boundaries with air or the 
third medium. It is straight-forward to take these boundary conditions into consideration. 
This procedure was discussed in detail in our previous paper [1]. 



V. EXPERIMENTAL POSSIBILITIES 

The variables for the S-matrix given in Sees. Ill and IV are determined by the optical 
parameters, namely the two phase-shifts and one reflection/transmission coefficient. The 
combinations of these three variables will determine the form of the S-matrix, which may 
take three different forms. 

We note first the that iV-dependence of the S-matrix comes from the form of A matrix 
or the RXR matrix of Eq.(40). If the optical parameters are in such a way that the A 
matrix takes the form of Eq.(43), the elements of the A N matrix of Eq.(51) are bounded and 
oscillating functions of N. If A takes the form of Eq.(45), the A N matrix becomes Eq.(52). 
The elements of this matrix are not bounded as N becomes large. Thus, in the real world, 
iV-layers can have two different types depending on the form of A. 

In addition, the optical layers can satisfy the condition that the expression of Eq.(47) be 
zero: 

sinh A — (sin a) cosh A = 0. (60) 

Then the RXR matrix takes the form of Eq.(48), and the N dependence is linear. This 
case can be tested as the optical parameters are varied from positive values of Eq.(47) to 
a positive value through zero. This condition does not depend on N. We have discussed a 
similar case in our previous paper [1]. 

In their recent paper [7], Baskal and Kim noted the same transition process for one- 
lens optics. They noted that the camera focusing mechanism corresponds to contraction of 
Wigner little groups. It is interesting to note that the same contraction mechanism exists 
in iV-layer optics. 



CONCLUDING REMARKS 

Based on Wigner's little groups, we have developed an algebraic method which allows us 
to study the cyclic properties of two-by-two S-matrices for multilayer optics. Starting from 
the single-layer S-matrix, it is possible to write the iV-layer matrix by multiplying one of 
the parameters by N. The iV-dependence is therefore transparent. 

This is possible because the core matrices of the Wigner's little groups have a sliderule 
property which allows us to perform multiplications by additions, as noted in Eq.(2). This 
property is an important element in computer designs. 
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As was noted in Ref. [7], the transition from Eq.(43) to Eq.(20) corresponds to camera 
focusing in one-lens optics. From the mathematical point of view, it corresponds to the 
contraction and expansion of the little groups. From the geometrical point of view, this 
corresponds to transformation from a circle to hyperbola. It is interesting to note that 
we can perform these operations also in multilayer optics. Indeed, as in the case of lens 
optics [7], multilayer optics can serve as an analogue computer for group contractions. 

The correspondence between the Lorentz group 0(3, 1) and SX(2,c), the group of two- 
by-two unimodular matrices, is well known. Since most of the matrices in ray optics are 
two-by-two, the Lorentz group is becoming the major language in this field. Ray optics is the 
backbone of future technology, and optical devices such as polarizers, lenses, interferometers, 
mutilayers, all speak the language of the Lorentz group. Thus, it is possible for the Lorentz 
group to play computational roles in future generations of computers. 

It is a prevailing view in physics, especially in optics, that group theory is only for 
studying symmetries and not useful for computational purposes. Indeed, we do not need 
group theory to carry out matrix multiplications given in this paper, and we started only with 
three matrices given in Eq.(17) and Eq.(18). However, are going through some important 
theorems in group theory while going through the simple matrix algebras given in this paper. 
We choose not to elaborate on this point. 
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